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Abstract 

The potential habitability of a terrestrial planet is usually defined by the possible existence of liquid water on its 
surface, since life as we know it needs liquid water at least during a part of its life cycle. The potential presence of liquid 
water on a planetary surface depends on many factors such as, most importantly, surface temperatures. The properties 
of the planetary atmosphere and its interaction with the radiative energy provided by the planet’s host star are thereby 
of decisive importance. 

In this study we investigate the influence of different main-sequence stars (F, G, and K-type stars) upon the climate 
of Earth-like extrasolar planets and their potential habitability by applying a state-of-the-art three-dimensional (3D) 
Earth climate model accounting for local and dynamical processes. The calculations have been performed for planets 
with Earth-like atmospheres at orbital distances (and corresponding orbital periods) where the total amount of energy 
received from the various host stars equals the solar constant. In contrast to previous 3D modeling studies, we include 
the effect of ozone radiative heating upon the vertical temperature structure of the atmospheres. The global orbital 
mean results obtained have been compared to those of a one-dimensional (ID) radiative convective climate model to 
investigate the approximation of global mean 3D results by those of ID models. 

The different stellar spectral energy distributions lead to different surface temperatures and due to ozone heating to very 
different vertical temperature structures. As previous ID studies we find higher surface temperatures for the Earth-like 
planet around the K-type star, and lower temperatures for the planet around the F-type star compared to an Earth-like 
planet around the Sun. However, this effect is more pronounced in the 3D model results than in the ID model because 
the 3D model accounts for feedback processes such as the ice-albedo and the water vapor feedback. Whether the ID 
model may approximate the global mean of the 3D model results strongly depends on the choice of the relative humidity 
profile in the ID model, which is used to determine the water vapor profile. Hence, possible changes in the hydrological 
cycle need to be accounted for when estimating the potential habitability of an extrasolar planet. 

Keywords: extrasolar planets, Earth-like atmosphere, climate, atmospheric dynamics, habitability 


1. Introduction 

As instrumental sensitivity increases the number of po¬ 
tentially rocky extrasolar planets detected has steadily in¬ 
creased since the 1990s. For these kind of objects the ques¬ 
tion of habitability, i.e. their potential to have liquid water 
on the planetary surface, is of special interest. Whether or 
not liquid water is in principle possible on the planetary 
surface depends, especially on surface temperatures and 
pressures, hence planetary climate. Therefore, from the 
atmospheric modeling point of view, a major objective is 
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to identify and understand the key processes determining 
the planetary climate. 

Knowledge of the nature of these small planets is rather 
limited, especially for those within the habitable zone. 
The classical habitable zone is given by the orbital dis¬ 
tances at which liquid water is possible on the surface of 
an Earth-like, rocky extrasolar planet assuming e.g. cer¬ 
tain atmospheric compositions and masses. The bound¬ 
aries of the habitabl e zone have been studied with ID mod ¬ 
els, e.g. Hart ( 1979 1. Kasting et al. (1993), Forget] ( 1998k 
von Paris et al.l ( 201 Bah . Pierrehumbert and Gaidos (2011), 


Kopparapu ct al. ( 2013k and recently also with 3D mod¬ 
els. e.g. lAbe et all ( 2011 ). Wolf and Toon (2014), Leconte 
et al. (2013b), Yang et al. (2014). For nearby transit- 
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ing planets a mean mass density may be determined via 
independent measurements of planetary radius, via the 
transit method, and planetary mass, via the radial ve - 
locity method, as e.g. for CoRoT-7b ( Leger et al.l . l2009h . 
Kepler-lOb (Batalha et ah, 2011, Dumusque et ah, 2014) or 
GJ1214b (Charbonneau et al., 2009). These kinds of inde¬ 
pendent measurements are unfortunately not yet available 
for potentially rocky planets within the habitable zone. 
However, for most of the planets and planetary candi¬ 
dates the orbital distance and the type of central star 
have been determined. This allows the estimation whether 
these planets lie within the so-called habitable zone. The 
influence of the stellar type upon planetary atmospheres 
of terrestrial planets has been studied ex t ensive l y wit h 
ID models, as e.g. in Easting et al.l (Il993l). Selsisl ( 2000 h . 
Segura et al.l ( 20031 20051) . iGrenfell et al. ( 2007 1. Kitzmann 


et al. (2010) Rauer et al. (2011), Kopparapu et al. (2013), 
Hedelt et al. ( 20131) . Rugheimer et al. (2013), showing that 


not only the total amount of stellar energy has a large 
impact on a planetary atmosphere but also its spectral 
energy distribution (SED). The strongly wavelength de¬ 
pendent absorption and scattering of stellar light by plan¬ 
etary atmospheres influences the surface temperature, the 
vertical temperature structure and also the atmospheric 
chemistry, which then may lead to different spectral ap¬ 
pearances. 3D modeling studies of terrestrial extrasolar 
planets have helped to gain insight into key processes and 
boundary conditions important for the climate of rocky 
extrasolar planets. These planets may be very different 
from Earth due to different rotation rates (Joshi et al., 
1997, Joshi 2003, Yang et al., 2013, 2014), different obliqui¬ 
ties (Williams and Pollard, 2003), eccentricities (Williams 
and Pollard, 2002 ), and water reservoirs (Abe et al., 2011, 
Leconte et al., 2013b). These 3D modeling studies all in¬ 
dicate that especially the hydrological cycle may have a 
large impact on planetary habitability. 

SEDs different from the Sun have been included in 3D at¬ 
mosphere stu dies, such asIW ordswort h et al.l (1201 ill. Lecont e 
et al (^Omh lShields et al.l (1201 3ll2014l) . lYa/ng et all (12013 . 
20141) . The influence of the different SEDs has however 


only been analyzed in detail bv IShiclds et al. ( 2013 . 2014 ). 
who investigate the influence of the stellar SED upon the 
ice-albedo effect. So far all 3D modeling studies including 
different stellar SEDs omit ozone in their atmospheres. 
Ozone (O 3 ) likely has a large impact on the vertical tem¬ 
perature structure of such atmospheres and thereby on 
the dynamical processes. In the search for O 3 as a biosig¬ 
nature, its effect on the temperature need s to b e take n 
into account . Pho tochemistry studies bv ISelsisI ( 2000l) . 
Segura et al.l ( 20031) and Grenfell et al. (2007) showed that 


an O 3 layer may be expected for planets around K- and 
F-type stars with Earth-like atmospheres. 

The impact of the SED upon the temperature structures 
of Earth-like atmospheres and the consequences for plane¬ 
tary climate have not yet been studied with a 3D climate 
model as previous studies neglect the influence of O 3 upon 
the temperature structures. It will be studied in detail 


here. 

This paper investigates the influence of the spectral 
stellar energy input upon Earth-like planetary atmospheres 
utilizing a state-of-the-art 3D General Circulation Model 
(GCM) of Earth. The response in temperatures, surface 
conditions and the hydrological cycle are analyzed and the 
global mean results are compared to those of a ID cloud- 
free radiative-convective climate model. 

Section [5] gives details about the models used, followed by 
a description of the scenarios studied in section D31 In the 
results section ([3]) we will show the influence of the SEDs 
upon temperatures m, the hydrological cycle (13.21) , and 
the surface conditions (13.31) . The sensitivity of the 3D 
model results to selected model parameters is discussed in 
section 1331 In section 1331 the 3D model results are com¬ 
pared to those of a ID climate model. The paper closes 
with a summary and conclusion (sec. [J). 

2. Computational details 

To analyze the impact of different stellar spectra upon 
the climate of Earth-like extrasolar planets and the impor¬ 
tance of processes such as the water vapor and albedo feed¬ 
back or atmospheric dynamics we use a 3D Earth climate 
model. The results are analyzed and the global orbital 
means are compared to those of a ID radiative-convective 
climate model. Details of the atmospheric models and the 
modeling scenarios are described in the following text. 

2.1. 3D atmospheric model for Earth-like planets 
For the 3D atmospheric model calculations of Earth-like 
exoplanets we use the EMAC (ECHAM/MESSy Atmo¬ 
spheric Chemistry) model (Jockcl et al., 2006), which has 
been developed for detailed investigations of the Earth’s 
climate. It uses the first version of the Modular Earth Sub¬ 
model System (MESSyl) to link multi-institutional com¬ 
puter codes. The core atmospheric model is the 5th gener¬ 
ation European Centre Hamburg GCM (ECHAM5, Roeck- 
ner et al., 2006). We apply EMAC (ECHAM5 version 
5.3.01, MESSy version 1.8) in T42L39-resolution, i.e. with 
a spherical truncation of T42 (corresponding to a quadratic 
Gaussian grid of approx. 2.8° by 2.8° in latitude and lon¬ 
gitude) and 39 hybrid pressure levels from the planetary 
surface up to 0.01 hPa, which corresponds to about 80 km 
for the Earth. Near the surface the grid is terrain fol¬ 
lowing whereas constant pressure levels are used in the 
upper atmosphere. The model setup features the key at¬ 
mospheric processes determining the planetary climate, 
i.e. radiative transfer, convection, the hydrological cycle 
including cloud processes. Atmospheric chemistry is ne¬ 
glected in this study, instead a fixed Earth-like atmosphere 
is assumed (see sec. [231). 

In EMAC the radiative transfer in the shortwave and 
the longwave regimes are treated separately. In the short¬ 
wave regime the radiative transfer of the stellar radiation 
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is calculated while in the longwave regime the transfer of 
the thermal radiation originating at the planetary surface 
and within the planetary atmosphere is treated. In the 
shortwave regime EMAC offers a high resolution scheme, 
FUBRAD (Freie Universitat Berlin high-resolution RA- 
Diation scheme), originally developed for solar variabil¬ 
ity studies (Nissen et ah, 2007), which operates at pres¬ 
sures lower than 70hPa, i.e. in the stratosphere and meso¬ 
sphere. In FUBRAD the absorption of stellar radiation by 
molecular oxygen (O 2 ) and O 3 is calculated in 49 bands 
ranging from 121 .4nm to 682.5nm. At higher pressures, 
i.e. at lower heights, the standard radiation scheme RAD4- 
ALL (Fouquart and Bonnel, 1980) is applied, which treats 
the radiative transfer in the UV and visible in one band 
ranging from 250 nm to 690 nm. The shortwave radiative 
transfer at near-infrared (NIR) wavelengths is calculated 
in three bands ranging from 690 nm to 4 pm over the en¬ 
tire vertical domain. Absorption by water (H 2 O), carbon 
dioxide (CO 2 ), O 3 , methane (CH 4 ), nitrous oxide (N 2 0), 
and carbon monoxide (CO), as well as Rayleigh scattering 
by air, scattering by aerosols, liquid and icy cloud particles 
are considered. The scheme uses the <5-Eddington approx¬ 
imation. 

For the calculation of the thermal radiation in the long¬ 
wave regime the Rapid Radiative Transfer Model (RRTM, 
Mlawer et ah, 1997) is used, which takes into account the 
thermal emission of the atmosphere and planetary surface, 
the absorption by radiative gases (H 2 0, CO 2 , O 3 , CH 4 , 
N 2 O, and chlorofluorocarbons (CFCs)), aerosols and cloud 
particles in 16 spectral bands (from 3.08 /im to 1000 pm) 
using the correlated-k approach. 

Convective transport of dry static energy, momentum, and 
moisture is calculated using the ECMWF (European Cen¬ 
tre of Mediu m-range Weathe r Forec asts) mass-flux convec¬ 
tion scheme ( Bechtold et all l2004h . The stratiform cloud 
coverage is calcu lated based on relative humidity following 
ISundovistl (1978). Cloud micro physics are parametrized 
following Lohmann and Roeckner (1996). The horizon¬ 
tal diffusion tendency is formul ated via a hyper -La placian 
based on the method bv lLaursen and Eliasenl ( 19891) . us¬ 
ing constant diffusion coefficients, which depend on the 
horizontal resolution. The vertical diffusion at the sur¬ 
face is obtained from a bulk transfer relation. Above the 
surface layer, eddy diffusion is assumed using diffusion co¬ 
efficients for moisture and heat, which are parametrized 
i n terms of turbulent kinetic energy and mixing lengths 
( Brinkop and Roeckneii IT 90 5>tl . 

At the lower boundary the atmospheric model is coupled 
to a mixed layer ocean model (Roeckner et al. 1995). It 
calculates the sea surface temperatures, sea ice coverage, 
and sea ice thickness for a mixed layer of 50 m from the 
net surface heat budget and a flux correction, the so-called 
q-flux. It accounts for the missing horizontal and vertical 
heat transport in the ocean and between the ocean and the 
atmosphere. The flux correction has been derived from a 
reference scenario, see section [2711 

Note that complex Earth climate models such as the one 


used here, rely on many complex parameterizations and 
parameters. Some of these parameters are used to ad¬ 
just the model in such a way that it properly reproduces 
present day and past climate states of the Earth. Espe¬ 
cially parameters concerning the cloud properties are often 
adjusted such as the asymmetry parameter of the water ice 
crystals or the number density of the condensation nuclei. 
This is one of the major reasons for differences in modeling 
results of Earth’s future climate, espe cially on seasonal and 
regional scales (see e.g. Flato et ah . 201311 . Furthermore, 
the vertical and horizontal resolutions could also have an 
influence affecting e.g. cloud covers and p recipitation pat¬ 
terns, as e.g. discussed for ECHAM5 by iRoeckner et al. 
(1200411 and f or th e Community Atmosphere Model 3 in 
Williamson ( 2008h . Also the choice of time steps and 


ti me scales may i nfluence the model results as shown by 
Williamsonl 02013 1. We usually apply a time step of 300 s 
but need to reduce it occasionally during spin-up to ensure 
small enough temperature tendencies. 

2.2. ID radiative-convective column model 
For comparison we compute vertical global mean atmo¬ 
spheric temperature and water profiles with a ID cloud- 
free radiative-convective column model, which ranges from 
the surface up to a height with a pressure of 6.6-10 ~ 2 hPa. 
The temperature profile is calculated from energy trans¬ 
port by radiative transfer and convective adjustment. The 
radiative transfer in the shortwave regime (from 237 nm to 
4.5 pm) is solved in 38 spect r al ban ds, using a <5-Eddington 
approximation ( Toon et'aO . 1 1989 1 and correlated-k expo¬ 
nential sums. In the longwave regime either RRTM is 
used, as in the 3D model calculations, or MRAC (Modified 
RRTM for App licatio n in CC> 2 -dominated atmospheres, 
von Paris et al. ( 201f)h ) for the comparison of different rel¬ 


ative humidity profiles, which includes only H 2 O and CO 2 
in 25 bands ranging from 1 to 500 pm. Whenever the lapse 
rate calculated via radiative equilibrium exceeds the adia¬ 
batic lapse rate, convective adjustment is performed to dry 
or moist adiabatic conditions. The water vapor profile is 
calculated from the temperature profile and a relative hu¬ 
midity (RH) parametrization by Manabe and Wetherald 
(1967). For the comparison of different relative humidi¬ 
ties (see section 15751) . additionally, a fully saturated atmo¬ 
sphere is assumed (RH=100%). A det ai led de s cripti on of 
the atmospheric mode l is given in I Hauer et al. ( 2011 ) and 
von Paris et al. ( 201flh and references therein. 


2.3. Modeling scenarios 

We investigate the influence of different main sequence 
stars upon the climate of Earth-like extrasolar planets. 
Planetary parameters are therefore chosen to resemble the 
present Earth, such as planetary radius, mass, hence grav¬ 
ity, obliquity, eccentricity, rotation rate, land-sea mask, 
and orography. Also an Earth-like atmospheric mass and 
chemical composition is assumed which is dominated by 
molecular nitrogen (N 2 ) and O 2 and includes spatially and 
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temporally uniform trace gas amounts of CO 2 (355ppm), 
CH 4 (1.64ppm) and N 2 0 (308ppb). For O 3 an annual 
mean of the zonal mean distribution given in Fortuin and 
Kelder (1998) is taken, see Fig. [TJ We assumed present 
day O 3 concentrations, as O 3 has a large impact on the 
temperature structure and thereby on the dynamical pro¬ 
cesses in our present day atmosphere. For the comparison 
of different relative humidities in section I5~5l an atmosphere 
composed only of N 2 , H 2 0 and C0 2 has been assumed. 

In the 3D model the amount of water in the atmo¬ 
sphere, i.e. water vapor, liquid water and water ice, is cal¬ 
culated considering various processes such as atmospheric 
transport, phase transitions within clouds and precipita¬ 
tion. In the ID model the amount of water vapor is de¬ 
termined by the assumption of a relative humidity profile 
and the calculated temperature profile. 
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Figure 1: Ozone concentrations (in ppmv) assumed in the 3D model 
calculations. 


In order to allow for ice free conditions in the 3D model 
calculations the background surface albedo of Antarctica 
and Greenland is set to 0.15 resembling the albedo of gran¬ 
ite. For the othe r glaci er free surface area the albedo map 
from lHagemann ( 2002f) is used. The surface albedo con¬ 
sidered by the radiative transfer is then calculated by the 
model taking the snow coverage of the land surfaces and 
canopy as well as the sea ice coverage and surface tem¬ 
peratures into account. However, no ice surfaces on land 
are calculated by the model, such as glaciers for instance. 
Therefore, the surface albedo for the reference scenario of 
the Earth around the Sun does not correspond to present 
Earth since snow on ice (e.g. glaciers) has a larger albedo 
than snow on granite. In the ID model calculations the 
surface albedo is set to 0.2 (or 0.22 for an atmosphere 
composed of N 2 , C0 2 and H 2 0 only, see sec. 1531) which 
is the value needed to obtain the mean surface temper¬ 


ature of the Earth (288 K) for a solar spect rum with a 
Total Solar Irradiance (TSI) of 1366Wm -2 ( Guevmard . 
200411 . to mimic the albedo effect of clouds. Note that 
the assumed land-sea mask, the surface albedo, as well as 
the water reservoir have an impact on the model results. 
For planets with a small water reservoir, the humidity of 
the atmosphere would be reduced, and different land-sea 
masks would also change e.g. the excitation of planetary 
waves, possibly leading to a different atmospheric circula¬ 
tion. The surface albedo is of great importance for the cal¬ 
culation of the radiative energy transport. In the present 
work, however we keep the number of parameters changed 
to a minimum as we are mainly interested in the impact 
of different stellar irradiations. 

For our central stars we use the Sun (G-type star) for 
the reference scenario, the F-type star er Bootis, and the 
K-type star e Eridani (see table [3}. The stellar spectra 
are depicted in Fig. [2] The orbital distances are chosen 
to yield a total energy input of 1366 Wm " 2 at the top of 
the atmosphere (TOA), which corresponds to the TSI of 
the present Sun. Since a Bootis is more luminous than 
the Sun this scaling of the total stellar energy flux cor¬ 
responds to a larger orbital distance and thereby longer 
orbital period. For e Eridani, which is less luminous, the 
opposite is the case. The orbital periods have been deter¬ 
mined using Kepler’s 3rd law and the stellar masses given 
in table |T| For these types of central stars and the rela¬ 
tively large orbital distances, one may reasonably assume 
that the planet s are not forced into synchronous rotation 
by tidal forces ( Kasting et ah . 1993h . therefore a rotation 
rate of the present Earth may be assumed. Table Q] summa¬ 
rizes the stellar properties as well as the orbital distances 
and periods of the planets. Details of the stellar spectra, 
which are a composite of satellite measurements in the UV 
and synt hetic model spectra at l arger wavelengths can be 


found in Ivitzmann et al. ( 201 Oh . The solar spectrum is 


taken from Guevmardl ( 2004 1. 


The mixed layer ocean model which is coupled to the 3D 
atmospheric model requires a flux correction, the so-called 
q-flux, to account for oceanic heat transport in the cal¬ 
culation of the sea surface temperatures. For Earth cli¬ 
mate calculations this heat flux is usually calculated from 
a reference scenario with prescribed sea surface tempera¬ 
tures (SSTs) and then applied to modeling scenarios with 
a small disturbance of the reference state. It therefore de¬ 
pends on the model and model setup such as horizontal res¬ 
olution. For large deviations from the reference state, a full 
atmosphere-ocean GCM is usually applied, e.g. for long¬ 
term Earth climate predictions. The utilization of a simple 
mixed-layer ocean model was preferred to a coupled atmo¬ 
sphere-ocean circulation model due to computational cost 
and the desire to limit the number of unknown boundary 
conditions. We therefore assume that the difference in the 
stellar spectral flux distribution to be a small disturbance 
of the reference state. The q-flux is calculated from net 
surface heat fluxes ( Fh ea t ), be. net radiative, sensible and 
latent heat flux, which have been calculated for a reference 
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scenario of the present Earth using prescribed climatolog¬ 
ical monthly mean sea surface temperatures ( Tsst ) from 
AMIP II (Atmospheric Model Intercomparison Project II, 
Taylor et ah, 2000) via: 


Q — Fheat Cm ^ , ( 1 ) 

with C m the heat capacity of the ocean. We make use 
of different q-flux corrections: 

• ql: the q-flux varies with every time step ( dt ). It 
uses monthly mean (mm) values of Fheat , which are 
calculated from a reference scenario, and dT Q ^ T , which 
varies with every time step and uses prescribed SSTs 

. This is the q-flux which would be applied for Earth 
climate calculations, as with this parametrization 
the prescribed SSTs of the reference scenario can be 
reproduced. 

For the calculation of planetary climates at different or¬ 
bital periods the above is however not useful, since it as¬ 
sumes a seasonality of a 365 days orbit. Therefore, differ¬ 
ent q-fluxes have been applied and their influence is dis¬ 
cussed in section E31 

• q2: the q-flux also varies with every time step but 
it uses the annual mean (am) of Fheat, and dT g^ T 
varies with every time step, 

• q3: the q-flux is the monthly mean (mm) of ql, 

• q4: the q-flux is the annual mean (am) of ql, and 

• q5: the q-flux is equal to zero. 

The q-fluxes q2, q3, and q4 represent only a small de¬ 
viation from the usually applied ql, mainly affecting the 
temporal consideration of the oceanic heat flux. Assuming 
the q5 q-flux completely ignores any heat transport. The 
q-fluxes q2-q4 test whther a small change in the oceanic 
heat flux can lead to large deviations of the climate stages, 
as expected for Earth climate simulations. The q-flux q5 
tests the robustness of our results. 

Note that for the results discussed in sections I3T1IO 
q2 is used for the planet around the K-type star and q4 is 
used for the planet around the F-type star. The difference 
in the annual global mean surface temperature owing to 
the different q-fluxes is however small as shown in section 

cm 


The seasons have different lengths depending on the 
orbital period. The length of the seasons, as e.g. north¬ 
ern hemispheric winter (NHW, used in sec. 13.11 and 13.211 . 
has been determined by the distribution of stellar insola¬ 
tion over latitudes and time. Hence, for each scenario the 
atmospheric temperature has been averaged over a time 
period for which the stellar insolation over latitudes cor¬ 
responds to e.g. December, January and February for the 



Figure 2: Stellar spectra (stellar energy flux Fg) used in this work. 
Blue: F-type star, Black: Sun, Red: K-type star. The spectra all 
yield a total energy input of 1366 Wm ' 2 at TOA, corresponding to 
the orbital distances given in tabled 


lame i: stellar ana ormtai parameters. 


Star 

stellar type 


a /au 

P /days 

a Bootis 

F2V 

1.194 a 

1.89 

868.64 

Sun 

G2V 

1 

1 

365.25 

e Eridani 

K2V 

0.82 b 

0.6 

184.00 


Bovaiian et al. ( 2012 1 b Butler et al. ( 2006lf 


Earth around the Sun. For the planet around the K-type 
star a NHW lasts 44 days and to 212 days for the planet 
around the F-type star. Furthermore, the results have 
been averaged over several orbits for which the atmosphere 
has reached a quasi-equilibrium state (after spin-up, when 
the influence of the initial state (e.g. stellar insolation) has 
ceased and surface temperature variations are caused by 
the seasonal cycle only): 6 orbits for the planet around 
the F-type star, 26 orbits for the planet around the Sun, 
and 18 orbits for the planet around the K-Star. For the 
planet around the F-type star this corresponds to more 
than 5000 days due to the longer orbital period (and only 
to about 3000 days for the planet around the K-type star). 
We limited the number of orbits included in the long term 
mean to save computing time. 

In section HOI the influence of different orbital periods, and 
hence different durations of seasons is discussed briefly. 
The orbital period has been varied because in principle 
different lengths of seasons may lead e.g. to different sea 
surface temperatures, because the ocean has a large ther¬ 
mal inertia. In addition, stellar parameters such as e.g. the 
stellar mass, which is used to calculate the orbital periods, 
also has uncertainties. This partly motivates the assump¬ 
tion of different various orbital periods. 
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3. Results and discussion 


To investigate the consequence of the different input spec¬ 
tra upon the climate we first discuss the results of the 3D 
model calculations for the Earth-like extrasolar planets or¬ 
biting different types of main-sequence stars in terms of 
surface temperatures, temperature structures, hydrologi¬ 
cal cycle and surface properties, and then compare these 
results to those of the ID model. 


3.1. Temperatures 

For the question of a planet’s habitability the tempera¬ 
tures at the surface are relevant. In the following first 
the response of the 2-meter-temperature (in the forthcom¬ 
ing text termed ’near-surface temperature’), i.e. the atmo¬ 
spheric temperature in the lowermost atmospheric layer, 
hence closest to the surface is discussed. Figure [3] depicts 
the orbital mean of the near-surface temperature for the 
Earth-like planets around different types of stars over lati¬ 
tudes and longitudes. Despite the fact that all three plan¬ 
etary scenarios receive the same amount of stellar energy, 
the near-surface temperatures are quite different. Temper¬ 
atures are lower for the planet around the F-type star and 
higher for the planet around the K-type star in compari¬ 
son to the Earth-like planet around the Sun. 


For the planet around the F-type star large regions of 
the planetary surface are ice free and habitable, with tem¬ 
peratures well above the freezing point of water (273.15 K) 
despite a global mean temperature of 273.6 K, see table 
[2] Temperatures below the freezing point of water are 
fou nd polewards of abo u t 40° latitude. For E arly Earth, 
e.g. iKunze et al. ( 2014h . Charnav et al.1 ( 201 .'ill . Wolf and 
Toon (2013) have also shown that global mean tempera¬ 
tures below the freezing point of water do not necessarily 
lead to freezing of the entire surface reservoir of liquid wa¬ 
ter, to so-called snowball states. Instead, regions of open 
water in the equatorial region may be found even for global 
mean temperatures as low as 250 K. This should be kept 
in mind when ID model results are used to evaluate the 
habitability of an extrasolar planet, Early Earth or Early 
Mars. 


For the planet around the K-type star the near surface 
temperature is everywhere higher than the maximum tem¬ 
perature obtained for the planet around the Sun, i.e. tem¬ 
peratures at polar latitudes for the planet around the K- 
type star are larger than equatorial temperatures for the 
planet around the Sun. The global orbital mean near sur¬ 
face temperatures (T 2 m ) for the scenarios are summarized 
in table |21 

These temperatures are the result of the interaction 
of various processes, such as absorption and scattering of 
stellar radiation by the atmosphere and surface, the green¬ 
house effect, and energy transport by convection and dy¬ 
namics. Our results confirm earlier studies (e.g. Kasting 



Figure 3: Orbital mean near surface temperatures (K) for Earth-like 
planets around the F-type star (upper panel), the Sun (middle panel) 
and the K-type star (lower panel). The thick black line corresponds 
to the melting temperature of water (273.15 K) 


et al., 1993, Segura et al., 2003, Kitzmann et al., 2010, 
Shields et al., 2013) which already discussed, that the 
spectral distribution of the stellar energy leads to differ¬ 
ent temperature responses of the atmosphere. The F-type 
star has its radiation maximum at shorter wavelengths and 
the K-type star at longer wavelengths than the Sun (see 
Fig. [2]). Therefore, radiation from the F-type star is more 
effectively scattered back to space via Rayleigh scatter¬ 
ing, while radiation by the K-type star is less effectively 
scattered and more effectively absorbed by water vapor, 
and carbon dioxide in the NIR. This leads to a higher 
planetary albedo for the planet around the F-type star 
and a lower planetary albedo for the planet around the 
K-type star (see table [2|. The planetary albedo has been 
determined from the outgoing longwave radiation assum- 
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Table 2: Global mean atmospheric and surface response for planets around different types of central stars. 


Star (Stellar type) 

T 2m (K) 

GHE (IC) 

Planetary albedo 
Water vapor column (■£&) 

Cloud cover (%) 

Cloud water column (-^) 

Cloud ice column |) 

Cloud GHE (K) 

Surface albedo 

Sea ice fraction relative to ocean (%) 
Snow depth (m) 


F3D 

Sun3D 

K3D 

a Boo (F2V) 

Sun (G2V) 

e Eri (K2V) 

273.6 

288.6 

334.9 

20.8 

39.2 

73.0 

0.38 

0.36 

0.23 

10.2 

30.1 

482.6 

0.75 

0.70 

0.67 

0.103 

0.106 

0.203 

0.037 

0.029 

0.028 

9.0 

8.9 

11.0 

0.21 

0.15 

0.10 

11.4 

3.9 

0.0 

0.0109 

0.0039 

0.0 


ing energy balance at the top of the atmosphere. Hence, 
the net incoming radiation (incident radiation - backscat- 
tered radiation) is thereby smaller for the planet around 
the F-type star and larger for the planet around the K-type 
star compared to a planet around the Sun, although the 
total top-of-atmosphcre (TOA) incident radiation flux is 
the same for all three cases. Furthermore more stellar ra¬ 
diation is absorbed by H 2 O and CO 2 for the planet around 
the K-type star, and less by the planet around the F-type 
star, respectively. This results in lower temperatures for 
the planet around the F-type star and in higher tempera¬ 
tures for the planet around the K-type star, as shown in 

Fig- 13 

This temperature response is intensified by positive cli¬ 
mate feedbacks, such as the water vapor feedback cycle or 
the ice-albedo feedback which strongly depend on local cir¬ 
cumstances, such as temperatures, but have an impact on 
the global climate state. 

For the planet around the F-type star the lower sur¬ 
face temperatures yield lower water vapor concentrations, 
which is discussed in section 13.21 see Fig. [8] The water 
vapor column is given in table [2] The lower water va¬ 
por concentrations lead to a decrease in the greenhouse 
effect. For the planet around the K-type star the opposite 
is the case, i.e. higher water vapor concentrations and a 
larger greenhouse effect. The greenhouse effect (GHE) in 
terms of temperature (given in table [2]) has been inferred 
from the difference of the temperature corresponding to 
the upwelling longwave radiation at the surface and the 
temperature corresponding to upwelling longwave radia¬ 
tion at TOA via the Stefan-Boltzmann law. The GHE for 
the planet around the Sun is a little higher than that of 
the Earth (about 33 K), which is probably caused by the 
assumption of the low surface albedo map leading to little 
higher temperatures and hence also higher water vapour 
concentrations. The strength of the greenhouse effect is 
illustrated in terms of the net longwave radiation ( Fi Wtnet ) 
close to the surface in Fig. [I] Fi w , net is the difference 


of the upward and downward directed longwave radiation 
(Flu i^net — Fi w ,down F/w; .?/p ) , hence negative values cor¬ 
respond to upward directed, positive values to downward 
directed radiation, respectively. Values close to zero, as 
for the planet around the K-type star, correspond to a 
large greenhouse effect, because a large amount of the up¬ 
ward directed radiation is absorbed and re-emitted by the 
atmosphere leading to less negative values. More nega¬ 
tive numbers close to the surface, as for the planet around 
the F-type star, correspond to a weaker greenhouse effect 
because the upward directed radiation is less effectively 
absorbed, and thereby the net longwave radiation is more 
negative. At larger heights the planet around the K-type 
star shows a more negative Fi w ^ net ^ due to the higher at¬ 
mospheric temperatures, and the planet around the F-type 
star less negative longwave radiation. 

The increase in temperature in the lower atmosphere 
for the planet around the K-type star is amplified by a cou¬ 
pling of the stellar radiation flux to the water vapor feed¬ 
back cycle. Due to the SED of the K-type star stronger 
radiative heating occurs in the lower atmosphere caused 
by the absorption of NIR radiation by water vapor, car¬ 
bon dioxide and clouds, as can be seen in the right panel 
of Fig. 0 This results in higher tropospheric and surface 
temperatures, leading to a higher amount of water vapor 
in the atmosphere by evaporation (see Fig. 0). This in¬ 
crease in water vapor enhances the greenhouse effect (see 
also Fig. 0 and the absorption of stellar radiation in the 
NIR, consequently yielding higher surface temperatures. 
This eventually leads to melting of all the ice and snow 
cover lowering the surface albedo (as discussed in section 
03 shown in Fig. m, which results in less reflection of 
stellar light at the planetary surface, hence more absorp¬ 
tion, causing even higher surface and tropospheric temper¬ 
atures. 

The different SEDs, however, not only result in a change 
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Figure 4: Global orbital mean net longwave radiation, F)iu n£t = 
Flu),down ~ Flw,up i f ()r the planet around the F-type star (blue), 
planet around the Sun (black), planet around the K-type star (red). 
The negative values indicate the upward direction. 




Figure 5: Global mean shortwave heating rates (for p<70 hPa) in the 
NIR (left) and due to O 3 (right) for the planet around the F-type 
star (blue), planet around the Sun (black), planet around the K-type 
star (red). 

of the surface conditions, but also affect the vertical tem¬ 
perature structure of the atmosphere. Figure [G] shows 
the zonal mean temperature structure for northern hemi¬ 
spheric winter (NHW, which corresponds to a December- 
January-February mean for the planet around the Sun) 
for all three scenarios from the surface up to a top pres¬ 
sure of 0.01 hPa. The zonal mean temperature structure 
for the southern hemispheric winter is not shown since the 
temperature differences between the hemispheres are small 
compared to the response to the different stellar spectral 
distributions. 

For the reference scenario of the planet around the Sun 
(middle panel of Fig. [ 6 ]) the typical vertical temperature 
structure of the Earth is visible with highest temperatures 
at the surface. Temperatures decline throughout the tro¬ 
posphere, reaching a minimum at the tropopause. Within 


the stratosphere temperatures increase reaching a maxi¬ 
mum in the stratopause and decrease again in the meso¬ 
sphere. The zonal temperature structure shows highest 
tropospheric temperatures in the equatorial region, and 
larger temperatures in the summer than in the winter 
hemisphere. Note that the resulting tropospheric tem¬ 
peratures at the summer south pole are higher than for 
the Earth, which is a result of the low background albedo 
map, see sec. m The tropopause is coldest in the equa¬ 
torial region and during polar night. At the stratopause 
highest temperatures occur during polar day due to the 
absorption of stellar radiation by O 3 . Similar tempera¬ 
tures at the stratopause are also found for polar night. 
These are the result of the meridional circulation, also 
called Brewer-Dobson-Circulation, which is indicated by 
the black contour lines in Fig. [G] (solid lines show clockwise, 
dashed lines anti-clockwise circulation, respectively). Air 
is lifted from the equatorial tropopause through the sum¬ 
mer stratosphere and mesosphere, and then transported to 
the winter hemisphere where it sinks. During ascent the 
air cools adiabatically, causing a cold summer mesosphere, 
and during descent the air heats adiabatically, causing the 
temperature maximum at the winter stratopause. 

For the planet around the F-type star (upper panel 
of Fig. ©, the overall temperature structure is similar to 
that of the Earth. However, tropospheric temperatures are 
lower than for the solar spectrum, whereas stratospheric 
temperatures are much higher, reaching maximum values 
of more than 350 K, i.e. larger than the surface tempera¬ 
ture. The high stratospheric temperatures in the summer 
hemisphere are caused by the high amount of stellar radi¬ 
ation at wavelengths where ozone absorbs, which is clearly 
visible in the ozone heating rates shown in the right panel 
of Fig. [5] The stratopause temperatures in the winter 
hemisphere are similar to those of the summer hemisphere. 
They result from the stratospheric meridional circulation, 
which is illustrated by the black contour lines in Fig. [G] 

For the planet around the K-type star (lower panel 
of Fig. [ 6 ]) the overall temperature structure differs con¬ 
siderably from that of the reference scenario around the 
Sun. There is no pronounced temperature increase within 
the stratosphere unlike the other two scenarios, because 
the stellar spectrum features much reduced radiation in 
the wavelength region where absorption by ozone occurs 
(at about 200-800 nm) compared to the Earth around the 
Sun, which results in less ozone heating (right panel of 
Fig. [5]). Furthermore, the pronounced cold tropopause in 
the equatorial region is missing. At lower pressures, how¬ 
ever, atmospheric temperatures become very low (below 
190 K) consistent with adiabatic cooling of air during as¬ 
cent. 

Due to the high surface temperatures the troposphere 
expands and shows nearly constant temperatures over all 
latitudes up to a pressure of 50 hPa. At pressures of about 
5 hPa features similar to the cold lower stratosphere of the 
Earth appear in both polar regions. Despite the fact that 















no pronounced stratospheric temperature maximum ex¬ 
ists in the summer hemisphere, a temperature maximum 
in the winter stratosphere is evident. This results from 
meridional circulation in the upper atmosphere from the 
equator to the poles (indicated by the black contour lines), 
which leads to a warming of the upper winter hemisphere 
due to adiabatic warming of descending air parcels. Note, 
that the meridional circulation patterns in the stratosphere 
are different for the planet around the K-type star com¬ 
pared to the planet around the Sun and the F-type star. 
The air rises from the equatorial troposphere through the 
equatorial stratosphere and mesosphere and sinks in both 
polar regions. This is in contrast to the Brewer-Dobson- 
Circulation where the air rises through the summer strato¬ 
sphere and mesosphere and sinks in the winter hemisphere. 

This change in the meridional circulation is a result 
of the different temperature structures which leads to a 
strong change in the zonal wind structure as shown in 
Fig. [3 For all scenarios we find westerly winds in the 
troposphere. These are caused by the temperature gra¬ 
dient from the equator to the pole. The warm air flows 
poleward and is deflected to the east due to the Coriolis 
force. For the planet around the F-type star the west¬ 
erly winds are weaker than for the planet around the Sun 
due to the smaller tropospheric temperature gradient. The 
stratospheres for the planets around the F-type star and 
the Sun show a westerly wind in the winter hemisphere 
and easterlies in the summer hemisphere, a result of the 
strong radiative heating of the summer stratosphere which 
leads to a reversed temperature gradient compared to the 
troposphere. Since the stratospheric temperature gradi¬ 
ents are larger for the planet around the F-type star the 
stratospheric zonal wind is stronger for this scenario. The 
distinct temperature increase above the summer pole is 
missing for the planet around the K-type star, thereby the 
latitudinal temperature gradient does not reverse at lower 
pressures and no easterly wind arises in the summer strato¬ 
sphere. Instead, westerly winds range over a large pressure 
region throughout the troposphere and stratosphere. The 
zonal wind plays an important role for the upward propa¬ 
gation of planetary waves, which force the meridional cir¬ 
culation. While planetary waves can propagate through 
westerlies, if the wind speed is not too large, they cannot 
pass through regions with easterly winds. For the planet 
around the Sun and the F-type star upward wave propa¬ 
gation is therefore damped in the summer hemisphere by 
the easterly wind leading to an asymmetric meridional cir¬ 
culation. The absence of a stratospheric easterly wind for 
the planet around the K-type star allows for the upward 
propagation of these waves, leading to a more symmetric 
meridional circulation of the upper atmosphere. A similar 
pattern as for the planet around the K-type star is found 
for the planet orbiting the Sun during springtime, where 
no strong easterly jet is present in the stratosphere. 

It is expected that the change in the temperature struc¬ 
ture would be even more pronounced when taking the 




Figure 6: Zonal mean atmospheric temperature (K) structure for 
northern hemispheric winter (NHW) for planets around different 
stars. Upper panel: planet around the F-type star, middle panel: 
planet around the Sun, lower panel: planet around the K-type star. 
The black contour lines show the meridional residual mass stream 
function in 10 9 kgm "" -, where solid lines represent clockwise, dotted 
lines anti-clockwise circulation. 


change in the atmospheric ozone concentrations due to dif¬ 
ferent stellar SEDs into account since ID modeling studies 
have shown that ozone concentrations are lower for Earth¬ 
like planets around K-type stars, which will lead to even 
lower stratospheric temperatures, and larger for planets 
around F-type stars resulting in even higher stratospheric 
temperatures (Selsis 2000, Segura et al., 2003, Grenfell et 
al., 2007). 
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Since previous 3D modeling studies of planets around 
other stars did not consider the heating by ozone, our 
study is the first to show the impact of the stellar spectral 
energy distribution upon the temperature structure and 
the atmospheric dynamics arising from ozone heating for 
planets around different host stars. While for the planet 
around the F-type star and the Sun this leads to very 
different temperature structures to those discussed in the 
literature, the temperature structure we find for the planet 
around the K-type star is similar to those f ound by other 
studies for atmospheres without ozone (e.g. IShields et all 
201,41201 4 . showing a decline of temperature with height. 
In the study presented here the lack of the stratospheric 
temperature increase leads to a change in the meridional 
circulation of the upper atmosphere for the planet around 
the K-type star. This will lead to a change in the trans¬ 
port of trace gas species such as ozone compared to the 
planet around the Sun, possibly influencing its total col¬ 
umn amount, as ozone is known to build up above the 
winter poles for the Earth. This possible change in the 
ozone transport, as well as the change in the temperature 
structure, will have an impact on e.g. the spectral appear¬ 
ance of such a planet. 


3.2. Hydrological cycle 

Water vapor plays an important role for the climate of 
a planet since it is an effective greenhouse gas. It can 
furthermore undergo phase changes, affecting the atmo¬ 
spheric temperature by the release or storage of latent 
heat. Water clouds scatter stellar radiation and contribute 
to the greenhouse effect. Furthermore, liquid water is a 
prerequisite for life as we know it. A planet with a liquid 
water reservoir on its surface is expected to also have water 
in its atmosphere which makes it a very important atmo¬ 
spheric constituent. Additionally the water vapor feedback 
can intensify climate responses: Higher surface tempera¬ 
tures cause a higher evaporation of water from the surface 
resulting in an increase of atmospheric water vapor and 
thereby a larger greenhouse effect and higher surface tem¬ 
peratures. However the overall effect is uncertain as the 
formation of clouds may alter the response as has been 
shown by Wolf and Toon (2014), Leconte et al. (2013b), 
Yang et al. (2013, 2014). The uncertain impact of clouds 
for other climates has been w idely discussed fo r climate 
change predictions, see e.g. IBouchcr et al.l (1201311 . 

Figure |S] shows the zonal mean distribution of the water 
vapor volume mixing ratio for the three scenarios in north¬ 
ern hemispheric winter. For the planet around the Sun the 
highest volume mixing ratios (vmrs) of around 0.01-0.02 
occur for near surface equatorial regions. From low to high 
latitudes the water vapor mixing ratio decreases steeply 
due to the decrease in temperature and hence evapora¬ 
tion. The summer hemisphere shows higher mixing ratios 
as expected. 

For the planet around the F-type star the water vapor 
mixing ratio is lower than for the reference scenario. Max¬ 
imum water vapor vmrs in the summer equatorial regions 
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Figure 7: Zonal mean wind (m/s) structure for northern hemispheric 
winter (NHW) for planets around different stars. Upper panel: 
planet around the F-type star, middle panel: planet around the Sun, 
lower panel: planet around the K-type star. Positive values indicate 
westerlies, while negative values show easterlies. 


are lower than 0.01. As for the planet around the Sun, 
water vapor concentrations steeply decrease towards the 
poles. Due to the lower water vapor amount in the at- 
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mosphere the greenhouse effect becomes less efficient, as 
indicated by the more negative near surface values of the 
net thermal radiation in Fig. [I] 


For the Earth-like planet around the K-type star much 
higher water vapor mixing ratios are obtained than for the 
planet around the Sun, reaching vmrs higher than 0.12 
in the equatorial region near the surface. Mixing ratios 
typical for the surface of the Earth occur at much higher 
altitudes at pressures of about 100-200 hPa. Here the con¬ 
centrations are nearly constant over all latitudes, which re¬ 
sults from a change in the hydrological cycle, which leads 
to recycling of water in the upper atmosphere. While 
for the Earth, cloud formation and precipitation lead to 
efficient replenishment of surface water, Figure [9] shows 
that for the planet around the K-type star a large part 
of the precipitation evaporates or melts before reaching 
the planetary surface. The high atmospheric temperatures 
lead to melting of all solid precipitation at pressures larger 
than 200 hPa. The amount of water recycled by this melt¬ 
ing and evaporation processes in the upper atmosphere 
constitutes a reservoir larger than the precipitation reach¬ 
ing the surface for the planet around the Sun. Nevertheless 
rainfall at the surface is intensified on the planet around 
the K-type star compared to the reference scenario. Thus, 
these high water vapor concentrations are the result of a 
change in the water vapor feedback cycle and the hydro- 
logical cycle. 

The high water vapor mixing ratios for the planet around 
the K-type star are also a result of the assumed land-sea 
mask, as about 70% of the planet assumed is covered with 
water. For a planetary scenario with a much smaller wa¬ 
ter reservoir, e.g. a planet without any ocean, the impact 
of the hydrological cycle upon surface tem perature s wo uld 
be much smaller as discussed e.g. by lAbe et ah ( 2011 ). 
since the greenhouse effect of water vapor would be much 
weaker in such a scenario. For a planetary scenario with 
a larger ocean reservoir, we expect a similar response of 
the hydrological cycle as the amount of water vapor in the 
atmosphere is related to the surface temperature via phase 
equilibrium. Due to a lower surface albedo for an ocean 
planet, surface temperatures would probably be higher. 

A comparison of the water vapor concentrations for 
the planet around t he K-type star with model results of 
Kasting et al.l ( 19931 ). who investigated the water loss limit 


as an inner boundary of the habitable zone shows that the 
stratospheric water vapor vmr (ca. 2xl0 -4 ) is still below 
their critical value of 3xl0” 3 . Thus, according to these re¬ 
sults a water reservoir of one Earth ocean should be stable 
over 4.5Gyrs when considering diffusion li mited escape of 
hydrogen to space for this scenario. iKasting and Po llack 
( 19831) however calculated atmospheric concentrations in 
the upper atmosphere of a Venus-like terrestrial planet, 
starting with stratospheric water vapor concentrations sim¬ 
ilar to those obtained here and found that an increase in 
the solar EUV flux by a factor of ten may result in strong 
hydrogen escape. Therefore, it is possible that even though 


the water vap or concentrations are below the critical value 
calculated bv lKasting et al.l ( 1993h a different stellar spec¬ 
tral energy distribution, especially with increased high- 
energy radiation may still lead to a water loss for such a 
scenario. 
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Figure 8: Zonal mean distribution of water vapor (vmr) for plan¬ 
ets around different stars in NHW. Upper panel: results for planet 
around the F-type star, middle panel: planet around the Sun, lower 
panel: planet around the K-type star. Note the linear pressure scale. 


In addition to water vapor, also clouds can have a 
strong impact on exop lanetary climate as shown e.g. by 
Kitzman n et al.l ( 20101) . We cannot judge whether the for- 
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Figure 9: Global orbital mean rain and snowfall for Earth-like planets 
around different stars. The left panel shows the rainfall profile and 
the right panel the vertical snowfall profile. The planet around the 
F-type star is depicted in blue, the planet around the Sun in black, 
and the planet around the K-type star in red. 


mation of clouds in Earth-like atmospheres of rocky extra¬ 
solar planets is captured by the cloud scheme included in 
our 3D model, since even for Earth climate predictions the 


impa ct of clouds varies over a wide range (IBo u chcr et al 


201 31). A first step towards the understanding of the cloud 


feedback for rocky extrasolar planets could be a model 
comparison for such scenarios. For the planet around the 
F-type star we find an increase in cloud cover and a de¬ 
crease in cloud cover for the planet around the K-type star 
compared to the planet around the Sun, see also table [2] 
We find the same behavior for the water ice column. For 
the cloud water column however we find an increase for 
the planet around the K-type star and a decrease for the 
planet around the F-type star, see tabled The cloud wa¬ 
ter column is about twice as high for the planet around 
the K-type star compared to the planet around the Sun. 
The greenhouse effect caused by clouds (Cloud GHE in 
table EJ is calculated from the difference in surface tem¬ 
perature and the effective temperature from the longwave 
radiation of the planets for cloudy and clear sky condi¬ 
tions. Our results suggest the cloud GHE is larger for 
the planet around the K and the F-type star than for the 
planet around the Sun. This indicates that the influence 
of clouds critically depends on cloud properties and at¬ 
mospheric temperatures, and cannot be easily estimated 
without sophisticated cloud modeling. The differences in 
cloud impact of different Earth climate models upon Earth 
climate pre dictions, but also studi e s of exoplanet s c enarios , 
such as e.g. Wolf and Toon ( 2014 ): Leconte et al. ( 2013b ). 
show that much more work is needed to understand the 
behavior and influence of clouds for the Earth and for ex¬ 
oplanetary climates. For the planet around the K-type 
star, the cloud GHE constitutes only about 15% of the 
total GHE, while for the planet around the F-type star 
it is about 43% and 22% for the planet around the Sun. 
Hence, for the warm scenario the cloud GHE seems to be 
of minor importance. 


3.3. Surface conditions 

The change in temperature is linked to a change in the sur¬ 
face albedo via the albedo feedback. Higher temperatures 
lead to melting of sea ice whereas lower temperatures will 
lead to its buildup. Figure [TU] shows the zonal orbital mean 
surface albedo for the three scenarios studied. At low and 
mid latitudes the albedo is the same for all scenarios. At 
mid to high latitudes the albedo differs depending on snow 
and ice coverage. 

For the planet around the K-type star no sea ice or snow 
is present, leading to a uniformly low surface albedo. 

The highest albedos occur for the planet around the F-type 
star at the north pole, which is a result of strong buildup 
of sea ice. In the southern hemisphere highest albedos oc¬ 
cur also for the planet around the F-type star at latitudes 
where sea ice is present. 

For the planet around the F-type star the ice-albedo 
feedback causes an increase in surface albedo due to buildup 
of sea ice down to latitudes of about 40°. Despite the fact 
that water vapor concentrations are low, i.e. the green¬ 
house effect of water is smaller than for the planet around 
the Sun, the increase in surface albedo is not large enough 
to cause a runaway glaciation. Note however, that de¬ 
pending on the stellar SED, the effect of the surface albedo 
may differ, since especially the albedo of snow is stron gly 
wavele ngth-dependent, as e.g. shown bv ljoshi and Haberlel 
( 20121) . In our study we do not account for the wavelength 
dependence of the surface albedo. Including it may in¬ 
crease the ice-albedo fe edback for the plane t around the 
F-type star as shown bv IShields et al.l 1 2013 ). 
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Figure 10: Zonal orbital mean surface albedo for planets around 
different stars. The planet around the F-type star is depicted in 
blue, the planet around the Sun in black, and the planet around the 
K-type star in red. 



3.4 . Sensitivity tests 

For an extrasolar planet, besides many other parameters, 
also the oceanic heat flux te r m is a -priori unknown. It 
was shown, by e.g. lYang et al.l ( 20131) , [Hu and Yan3 ( 20131) 
and ICullum et al.l ()2014h , that the oceanic heat flux may 
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strongly influence atmospheric GCM results and surface 
temperatures. Therefore, we test the sensitivity of our 
model results to a change in the q-flux term. Additionally, 
we vary the orbital period to estimate its influence on the 
model results, since different lengths in season may lead 
to different results due to the large thermal inertia of the 
ocean. 

The scenarios and their resulting global orbital mean tem¬ 
perature are summarized in tabic [3] Overall, the difference 
in temperature due to different orbital periods is small, 
about 1 K for the planet around the F-type star (compar¬ 
ing scenarios F3D and F3Dq4, which both use the q4 q-flux 
but different orbital periods). However, the seasonality 
changes, hence for longer orbital periods, longer temper¬ 
ature maxima and minima are obtained, see Fig. Ilf I For 
the planet around the F-type star and the Sun the season¬ 
ality is more pronounced than for the planet around the 
K-type star. While for the planet around the K-type star 
the atmosphere is dominated by the large amount of water 
vapor, for the planet around the F-type star and the Sun 
the build-up and melting of sea ice are important factors 
for the seasonal temperature variability. Note that for the 
planet around the F-type star we conducted most of the 
sensitivity tests with an orbital period of 450 days, which 
was computationally less expensive than using the physi¬ 
cally correct orbital period of 868 days used for the results 
in the previous sections . 

The influence of varying the oceanic heat flux (q-flux), 
used in the mixed layer ocean model to calculate the sea 
surface temperatures and sea ice, is a little larger, lead¬ 
ing to a temperature difference of up to 2 K for moderate 
changes in the q-flux term. This term partly represents 
the oceanic circulation and heat transport. We applied 
five different q-flux terms ql-q5, as described in sec. 12.31 
Figure[TT]suggests that moderately varying the q-flux term 
(shown are q2-q4) does not change the global temperature 
response. A seasonal variation of the temperature can be 
identified for the planet around the Sun and the F-type 
star. For the planet around the K-type star the the varia¬ 
tion of the temperature is less pronounced. A larger vari¬ 
ability at smaller timescales is visible which is less periodic. 
This is caused by the large amount of water compounds 
in the atmosphere which are highly variable on timescales 
smaller than the seasons. 

While on the global scale the near surface temperatures 
deviate by up to 10 K for different q-fluxes and orbital pe¬ 
riods (except for the planet around the F-type star with 
a q-flux of 0 ), on the regional scale the temperature dif¬ 
ferences are larger. For the planet around the F-type star 
we find a temperature difference of about 30-50 K at the 
summer pole and of about 20 K for the winter pole due 
to the difference in orbital period, with higher tempera¬ 
tures in summer and lower temperatures in winter for the 
longer orbital period (not shown). For the planet around 
the Sun the different q-fluxes ql-q3 lead to a small change 
in the global annual mean near surface temperatures from 


an exoplanet perspective. For Earth climate calculations 
these differences of about 2K are however large. On the 
local scale we find e.g. larger near surface temperatures (of 
about 5 K) above the equatorial ocean and the southern 
pole for the q-flux q3. The zonal seasonal mean tempera¬ 
ture structures are not influenced strongly by the change 
of the orbital periods or q-fluxes either, from an exoplanet 
perspective. Compared with the scenarios with an orbital 
period of 365 days and the q-flux ql we find a very simi¬ 
lar temperature structure, hydrological cycle, and a simi¬ 
lar dynamical behavior, and surface response compared to 
those presented in the previous sections. It was not antici¬ 
pated that the change in the orbital period and the oceanic 
heat flux would have such a small effect on the mean 3D 
model results. 


Only the assumption of an extreme and unrealistic q- 
flux of zero (q5) leads to a strong change in surface condi¬ 
tions for the planet around the F-type star, which under¬ 
goes global glaciation. Hence, for this scenario, an oceanic 
heat flux is required to prevent the planet from freezing 
completely. For the planet around the K-type star the 
orbital mean temperature only decreases by about 5 K for 
the q5 q-flux, since the greenhouse effect of the atmosphere 
do minates t he heat bu dget. 

In Shields et alJ ( 2013h a similar scenario has been studied: 
an Earth-like planet around an F-type star with an oceanic 
heat flux of zero. They discuss that their planet is close 
to global glaciation for a total amount of incoming stellar 
radiation as Earth receives from the Sun. Hence, although 
they include a wavelength-dependent ice albedo in their 
study, which should enhance the response by increased 
scattering of stellar light, in their model calculations the 
planet does not undergo global glaciation. Their model re¬ 
sults differ from ours because they disregard ozone in their 
calculations and assume a surface completely covered with 
water which lowers the background surface albedo. Fur¬ 
thermore, the albedo they assume in the visible wavelength 
regime ( 0.8 cold dry snow) is about equal to the one we 
assume in our model calculations for the entire shortwave 
regime (0.8 for snow on ice). Their albedo in the NIR 
regime is however lower (0.68 for snow on ice). Addition¬ 
ally, from their paper it is unclear whether they also use 
different amounts of greenhouse gases, such as CO 2 and 
CH 4 , which may also yield higher surface temperatures. 


3.5. Comparison of the 3D to the ID model results 
To evaluate the habitability of rocky extrasolar planets 
ID atmospheric models are often utilized since the ID as¬ 
sumption keeps the number of boundary conditions small. 
Furthermore, in future, averaged quantities of these dis¬ 
tant worlds will be the first to be retrieved. However, it 
has been suggested that for rocky extrasolar planets close 
to the inner edge of the habitable zone 3D phenomena 
may have to be taken into account for certain scenarios 
of e.g. tidally-locked extrasolar planets or for planets with 
small water reservoirs (Abe et ah, 2011, Leconte et ah, 


13 








Table 3: Sensitivity test scenarios and global orbital mean near-surface temperatures. 


Scenario 

Star (stellar type) 

orbital period (d) 

q-flux term (variability) 

T 2m (K) 

F3D 

a Boo (F2V) 

868.64 

q4 (am) 

273.6 

F3DqO 

o Boo (F2V) 

450.0 

q5 (=0) 

234.0 

F3Dq4 

o Boo (F2V) 

450.0 

q4 (am) 

274.7 

F3Dq3 

(7 Boo (F2V) 

450.0 

q3 (mm) 

275.5 

F3Dq2 

cr Boo (F2V) 

450.0 

q2 

275.2 

F3Dql 

cr Boo (F2V) 

365.25 

ql 

275.0 

Sun3D 

Sun (G2V) 

365.25 

q2 

288.6 

Sun3Dq3 

Sun (G2V) 

365.25 

q3 (mm) 

289.7 

Sun3Dql 

Sun (G2V) 

365.25 

ql 

287.7 

K3D 

e Eri (K2V) 

184.00 

q2 

334.9 

K3DqO 

e Eri (K2V) 

184.00 

q5 (=0) 

329.2 

K3Dq2 

e Eri (K2V) 

184.00 

q3 (mm) 

336.2 

K3Dql 

e Eri (K2V) 

365.25 

ql 

334.5 



Figure 11: Influence of different q-fluxes and orbital periods on near¬ 
surface temperature for planets around different stars. Results for 
the planet around the F-type star in blue, for solar radiation in black, 
and for the planet around the K-type star in red. q-flux corrections: 
ql dot-dot dashed, q2 solid, q3 long dashed, q4 dot-dashed for the 
F3Dq4 scenario with a 450 days orbit, q4 long dashed-short dashed 
for the F3D scenario with a 868 days orbit, q5 dotted. 


2013b, Yang et al., 2014). 

Here we test whether 3D phenomena need to be taken 
into account for habitable Earth-like extrasolar planets 
which are not tidally locked by comparing the global or¬ 
bital mean results of the 3D model with those of a steady- 
state radiative-convective model (see section l2T2l) . The up¬ 
per panel in Fig. [T2] shows the temperature-pressure pro¬ 
files for the three scenarios as computed with the 3D and 
the ID models. 

For the planets around the Sun and around the F-type 
star, the temperature profiles compare well for pressures 
between 1000 hPa and «1 hPa. Differences at larger height 
may be explained by the different model regimes. In par¬ 


ticular, the 3D model extends to 0.01 hPa while the ID 
model only ranges to 0.066 hPa. Also in the ID model, 
e.g. the absorption of radiation at 121.5 nm (Lyman a) by 
molecular oxygen is not considered leading to lower tem¬ 
peratures in the upper model atmosphere. For the planet 
around the K-type star a larger difference between the 
global orbital mean profile of the 3D model calculations 
and the ID model can be found. This mainly results from 
the large difference in the calculated water vapor profiles, 
which are shown in the lower panel of Fig. [12] In the 
ID model the water vapor profile is determined from an 
assumed vertical profile of relative humidity. For the ID 
model calculations shown in Fig. Q2] (FID, SunlD, KID) 
a relative h u midity profi l e as m easured for present Earth 
( Manabe and Wetheraldl . [ l967 ) is assumed for the calcu¬ 
lation of the water vapor profile, which starts with a rel¬ 
ative humidity of 77% at the surface and then decreases 
throughout the troposphere. This results in the steep de¬ 
cline with height of water vapor in the atmospheres as 
calculated by the ID model. However, in the 3D model 
calculations, as shown in section EZU the distribution of 
water vapor in the atmosphere differs due to a change in 
the hydrological cycle. Therefore, the water vapor profile 
and thus also the temperature profiles as well as the global 
mean surface values (Tab. [4]) calculated by the ID and the 
3D model disagree for the planet around the K-type star. 


To test whether the 3D global orbital mean profiles can 
nevertheless be approximated by the ID model, we addi¬ 
tionally compared the 3D model results for the K-type star 
with ID mo del profi l es calc ulated with a relative humidity 
of 100% by IStracke] ( 2012h (K1DRH100). This study uti¬ 
lized the same atmospheric model however with e.g. an 
improved thermal radiative transfer sch eme (ap pli cable 


to larg er temperature and pressure range, Ivon Paris et al 


( 2010 1). and applied it to an Earth-like planet around dif¬ 
ferent types of stars however neglecting the influence of 
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(K1DRH100) the tropospheric profiles of the 3D model re¬ 
sults for the pla net aroun d the K-type star and ID model 
calculation from Strackel ( 2012 ) agree well. Surface tem¬ 
peratures and water vapor volume mixing ratios also com¬ 
pare better for this assumption, see table 2] 



Figure 13: Comparison of the global orbital mean tropospheric tem¬ 
perature profile of the 3D model results for the planet around the 
K-type star (red, solid) with ID model results llStrackd . 1201 2) for 
the corresponding scenario with diffe rent RH profiles: orange, dotte d 
with a relative humidity profile from lManabc and Wethcraldl if1907 ). 
green, short-dashed with a relative humidity of 100 % 


Figure 12: Mean temperature (top) and water vapor (bottom) pro¬ 
files for Earth-like planets around different types of stars. Global 
orbital mean results from the 3D model in solid, profiles from the 
ID model in dashed lines. Results the planet around the F-type star 
in blue, for solar radiation in black, and for the planet around the 
K-type star in red. 


Table 4: ID and 3D global mean temperature (T) and water 
vapor vmr (H 2 O) at 1000 hPa for the Earth-like planets around 
different stars as calculated in this study, iSfarlD/3D, and in 
Stracke (2012) with diffe r ent re lative humidity profiles: RHMW from 
iManabe and Wetheraldl il967l ) RH100 with a relative humidity of 
100 %. 


Scenario 

T(K) 

H 2 0 (%) 

F3D 

274.3 

0.6 

FID 

280.1 

0.8 

Sun3D 

289.0 

1.4 

SunlD 

287.7 

1.3 

K3D 

333.4 

15.6 

KID 

292.9 

1.7 

K1DRHMW 

292.0 

1.7 

K1DRH100 

329.0 

14.1 


oxygen, ozone, methane, and nitrous oxide. For these 
model calculations (K1DRHMW, K1DRH100) a surface 
albedo of 0.22 was assumed. As can be inferred from 
Fig- GUI for the assumption of a fully saturated atmosphere 


This implies that the ID model can indeed approximate 
the global orbital mean climate states calculated with the 
3D model, when assuming an appropriate relative humid¬ 
ity profile. Note, that assuming a relative humidity of 
100 % in the ID model increases the surface temperature 
of the planet around the Sun to 305 K (for an atmosphere 
with N 2 , CO 2 , and H 2 O). Thus the agreement of the 3D 
and ID model results strongly depends on the choice of 
the relative humidity profile used for the water vapor pro¬ 
file calculation. Which profile to choose is however a-priori 
unknown. 

Note that the realization of a global mean relative hu¬ 
midity of 100% seems unlikely as atmosp heric dynamics 


will lead to drying of air during ascent (jLeconte et al 


2013a ). by condensation and cloud formation for most sce¬ 


narios. In the 3D model calculations presented in sec. 13.11 
IQ we find a global mean surface relative humidity of 
only 74% for the planet around the K-type star (71% 
for the planet around the Sun and 69% for the planet 
around the F-type star), lower than the relative humid¬ 
ity of 100% assumed in the ID model to find a good 
agreement between the results. However, the decrease 
in relative humidity with height is weaker for the planet 
around the K-type star than for the planet around the 
Sun and the F-type star in the 3D model results. For 
these cool er scenarios the decreas e is si milar to the RH 
profile bv IManabe and Wetheraldl (1967) assumed in the 
ID model. Including the global mean relative humidity 
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profile from the 3D model calculations in the ID model 
results in an increased surface temperature which reaches 
315 K. We find a similar good agreement between the ID 
and 3D model results as with a constant RH of 100% in 
the ID model, when assuming a fixed water vapor profile 
corresponding to the 3D global mean water vapor profile 
and a global mean surface albedo of 0.1 in the ID model 
calculations for the planet around the K-type star, yield¬ 
ing a surface temperature of 328 K. 


From our 3D model calculations which include the change 
in the hydrological cycle, the water vapor feedback, and 
the ice-albedo feedback, we deduce a maximum difference 
in the global orbital mean near-surface temperature of 
61K, for the planet around the F-type star and the planet 
around the K-type star (see table Q| when including an 
Earth-like oceanic heat transport. Including also the q-flux 
sensitivity test scenarios (see section 13.411 the maximum 
difference is 102 K for the same amount of total energy 
incident at top of the atmosphere. The largest difference 
in near surface temperatures obtained by the ID model 
calculations in this work is 49 K, for the planet around 
the F-type star with a relativ e humidity profile following 
Manabe and Wetherald ( 1967h and the planet around the 


K-type star with a relative humidity of 100%. The maxi¬ 
mum temperature difference obtained by the ID model cal¬ 
culations is smaller because a change in the surface albedo 
is not captured by the ID model calculations. Note that 
the 3D model results depend on complex parameteriza- 
tions which are adjusted to reproduce present and past 
Earth climates. Therefore, the results may differ in ab¬ 
solute numbers for varying the parameter sets and model 
codes. Our results should hence be considered to be of 
more of a qualitative than of quantitative nature. 

In ID model studies, similar to that presented here, the 
surface albedo is often tuned to fit the mean temperature 
of the Earth (288 K), and therefore not only includes the 
real surface albedo, but additionally mimics the reflectiv¬ 
ity of clouds within the atmosphere. The surface albedo in 
the 3D model is calculated and differs for the scenarios in¬ 
vestigated here as discussed in section HT~3l As long as the 
planets stay habitable the difference in surface albedo is 
however relatively small (Aalbedo ~ 0.1). When including 
the sensitivity scenario where the planet around the F-type 
star undergoes global glaciation the difference in albedo is 
large (Aalbedo ~ 0.5). The impact of clo uds on the plan¬ 


etary albedo can be very diverse (see e.g. iKitzmann et a l 


etary 

120101 . and therefore difficult to capture. 

Which surface albedo is realized for a certain plan¬ 
etary scenario needs to be calculated interactively by a 
coupled 3D atmosphere-ocean general circulation model 
(AOGCM), since the resulting surface albedo strongly de¬ 
pends on the oceanic heat transport as illustrated in sec¬ 
tion [XT] Therefore atmospheric GCMs coupled to a mixed 
layer ocean, as utilized in this work, may only give a 
first estimate of the surface albedo for a certain planetary 
scenario. It would preferably be calculated with coupled 


atmosphere-ocean models. Such calculations will however 
introduce an even larger amount of parameters which are 
unknown for extrasolar planets, and are more expensive 
in terms of computing time. Their impact should never¬ 
theless be explored. It should be noted that the planetary 
surface albedo is important to deduce the climate of plan¬ 
ets with relatively thin atmospheres, as for these planets 
it has large impact on the energy budget. For planets 
with thicker atmospheres e.g. at the inner edge of the 
habitable zone, w here one may expect water dominated 
atmospheres (e.g. Strackd . 2012h or at the outer edge of 
the habitable zone where one may expect thick C02-do- 
minated atmospheres its impact is much smal ler because 


it is masked by the thick atmospheres (Ivon Paris et al 
2013blIShields et allliobill . 


4. Summary and conclusion 

The results of 3D GCM climate modeling of Earth-like 
planets around F, G, and K-type stars have been dis¬ 
cussed in this paper. It has been shown that different 
stellar spectral energy distributions may lead to very dif¬ 
ferent climate states of an Earth-like planet for the same 
total amount of stellar energy incident at the top of the 
atmosphere. This results from the wavelength-dependent 
absorption and scattering properties of the atmosphere, 
and an amplification of the climate response by positive 
climate feedback cycles, namely the water vapor and the 
ice-albedo feedback. 

For the planet around the F-type star we find an enhanced 
heating of the stratosphere by ozone due to the change 
in the stellar energy distribution compared to the planet 
around the Sun. For the planet around the K-type star we 
find a strong change in the vertical temperature structure 
with no temperature increase in the stratosphere, caused 
by decreased heating by ozone absorption, and a change in 
the hydrological cycle due to the high tropospheric tem¬ 
peratures. For nearly all scenarios studied, the Earth-like 
planets result in habitable surface conditions. Only when 
neglecting the (parametrized) oceanic heat transport for 
the planet around the F-type star, which is rather un¬ 
physical, we find uninhabitable surface conditions. For 
this scenario the planet undergoes global glaciation. A 
moderate change in the oceanic heat flux term and in the 
orbital period do not show a strong impact on the mean 
planetary surface climate. 

Comparing the global orbital mean of the 3D model re¬ 
sults to those of a cloud-free ID radiative-convective col¬ 
umn model showed that the temperature response may 
be approximated by the ID model. The agreement of the 
results, however, crucially depends on the choice of the rel¬ 
ative humidity profile utilized to calculate the water vapor 
profile. 
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